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O ■ 

^vq , We present a practical way of smoothing discrete breakup S-matrix elements calculated 
by the continuum-discretized coupled-channel method (CDCC). This method makes the 

' smoothing procedure much easier. The reliability of the smoothing method is confirmed for 

1} ' the three-body breakup reactions, 58 Ni(d,pn) at 80 MeV and 12 C( 6 He, 4 He 2 n) at 229.8 MeV. 
U-i ' 

O §1. Introduction 

Secondary beam experiments are opening up new physics on unstable nuclei 
away from the stability line in the nuclear chart. Such unstable nuclei have exotic 
properties, e.g., the halo structure. x )~ 3 ) In scattering of such a weakly bound pro- 
jectile, it easily breaks up into its constituents, and hence the reaction should be 
described as at least a three-body scattering. 

One of the most reliable methods for treating the projectile breakup processes 
in the wide range of incident energy is the method of continuum-discretized coupled 
channels (CDCC). 4 )' 5 ) In CDCC, the scattering wave function of the total system 
is expanded with the complete set of bound and continuum states of the projectile. 
\f} [ The continuum states are classified by linear and angular momenta, k and £, and 

truncated at upper bounds, /c max and £ max , respectively. This truncation is based 
{NJ ■ on the assumption that the breakup processes take place in the truncated space. 

The /c-continuum is divided into small bins and the continuum states in each bin 
are averaged into a single state. This discretization procedure is called the aver- 
age discretization method. The space spanned by the bound states and the finite 
number of discretized continuum states are called the modelspace. The 5-matrix 
elements calculated with CDCC converge as the modelspace is extended. 6 )' 7 ) The 
converged CDCC solution is the unperturbed solution of the distorted Faddeev equa- 
tions, and corrections to the solution are negligible within the spatial region in which 
the breakup processes take place. 8 )' 9 ) 

The average discretiztation method has widely been used so far, but its ap- 
plication is limited in the three-body breakup processes in scattering of two-body 
projectiles. The average method is not feasible for the four-body breakup processes 
in scattering of three-body projectiles, since it requires exact three-body continuum 
states of the projectile that are quite difficult to obtain. This problem can be solved 
by using the pseudostate discretization method 10 ) -16 ) in which the continuum states 
{■ifj(k)} are replaced by pseudostates {'Pi} obtained by diagonalizing the internal 
Hamiltonian of the projectile in a space spanned by L 2 -type basis functions. One 
can adopt the Gaussian 10 )' 11 ) or the transformed harmonic oscillator (THO) 12 ) ba- 
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sis as the L 2 -type basis functions. The validity of the pseudostate discretization 
method was confirmed in scattering of two-body projectile by the good agreement 
between the CDCC solutions obtained by the pseudostate discretization and the 
average methods. 10 )' n )' 13 )~ 15 ) 

For four-body reaction systems, CDCC calculations with the pseudostate dis- 
cretization method based on Gaussian 13 -* -15 ) or THO 16 ) basis functions were success- 
ful in describing the elastic scattering at not only high energies but also low energies 
near the Coulomb barrier. Thus, the back-coupling effects of four-body breakup 
processes on the elastic scattering, i.e., virtual four-body breakup processes, are well 
described by CDCC based on the pseudostate discretization method. 

The ^-matrix elements calculated with CDCC, Si, are discrete in k, while the 
exact ones S(k) are continuous. Thus, one needs a way of smoothing Si in order 
to analyze real breakup processes, i.e. breakup reactions themselves. In principle, 
this can be possible by calculating the overlap ^(k) = (0(fc)[(Pj) between the exact 
continuum states ip(k) and the pseudostates see H2.'2\ for the detail. In practice, 
this smoothing factors J~i{k) can easily be obtained in three-body breakup reactions 
of two-body projectiles, but it is quite hard in four-body ones of three-body projec- 
tiles. Very recently, we have proposed a smoothing method applicable to four-body 
breakup reactions, 17 ) but it still requires heavy numerical calculations and hence 
not so practical. Thus, it is highly expected that an accurate and practical method 
of smoothing Si is proposed. 

The complex scaling method (CSM) 18 ) is a powerful tool of solving many-body 
resonance and weakly bound states. Recently, CSM was extended so that it can 
treat the electromagnetic transition of the core+nucleon+nucleon system such as 
6 He and n Li from its ground state to the continuum ones. 19 )' 20 ) As a merit of 
CSM, the calculated transition strength is obtained as a continuum spectrum. This 
implies that an accurate and practical method of smoothing Si can be constructed 
with CSM. 

In this Letter, we present a simple formula that makes it possible to smooth Si 
accurately, using CSM. This formula is useful for not only three-body but also four- 
body breakup reactions. Validity of the formula is tested in the three-body breakup 
reactions, 58 Ni(d,pn) at 80 MeV and 12 C( 6 He, 4 He 2 n) at 229.8 MeV, since in such 
three-body reactions the "exact" S'-matrix elements are obtainable by calculating 
jFj(fc) exactly. This test is inevitable to proceed to CDCC calculations for four-body 
breakup processes in future. 

This Letter is constructed as follows. In Sec. [2j we first recapitulate CDCC 
and the pseudostate discretization method, and present an accurate and practical 
smoothing formula of Si. Test calculations are done in Sec. [3l Section 2] gives a 
summary. 
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§2. Formulation 

2.1. CDCC and the pseudostate discretization 

For simplicity, we consider the three-body (A+b+c) system, where particles 
A, b and c are assumed to be structureless. It is straightforward to extend the 
present formulation to the four-body system in which the projectile consists of three 
constituents. The total Hamiltonian of the A+b+c system is 

H = T R + U + H B , (1) 
U = U b (R b ) + U C (R C ) + V Con \R), (2) 
H B = T r + V bc , (3) 

where R is the relative coordinate between the center-of-mass (cm.) of the projectile 
B=b+c and A. Coordinate R x (x = b and c) denotes the relative coordinate between 
x and A, and T% (£ = R and r) is the kinetic-energy operator associated with £. 
Interaction Vt> c is the potential between b and c, and U x is the nuclear part of 
the optical potential between x and A. Meanwhile the Coulomb interaction y Coul 
between B and A is treated approximately as a function of R only, i.e., we neglect 
Coulomb breakup processes. 

In the pseudostate discretization method, 10 )' n )' 13 )~ 15 ) the pseudostates 
are obtained by diagonalizing the internal Hamiltonian H B in a space spanned by 
the L 2 -type Gaussian basis functions: 21 ) 

(&i\H B \&i>) = Swii, (4) 

where £j is the eigenenergy of d>j. The subscript i denotes a set of quantum numbers, 
i.e. the energy index n of the pseudostates, the angular momentum I between b 
and c, and its projection on the z-axis m. The eigenstates of H B with negative and 
positive energies correspond to the bound state(s) and the pseudostates, respectively, 
and the latter are regarded as discretized continuum states. 

The basic assumption of CDCC is that the breakup processes take place in a 
model space 14 )' 15 ) 

r = J2\$i)(H (5) 

i 

Validity of this assumption is justified by the fact that the calculated elastic and 
total breakup cross sections of the three- and four-body scatterings converge as the 
modelspace is extended. 10 )' n )' 13 )~ 15 ) Thus, one may regard {^} as a complete 
set in describing the reaction processes considered. We henceforth call {^i} the 
approximate complete set in this meaning. 

We solve the three-body Schrodinger equation in the modelspace V, 

V[H - E tot ]V\* CDCC ) = 0, (6) 



where E to t is the total energy of the system. The total wave function i^ CDCC is 
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expanded by the approximate complete set {^j}: 



i^CDCC\ 



(7) 



where \$i,Xi) = \<&i) ® \xi) and i = denotes the elastic channel and others (i ^ 0) 
the breakup channels. The expansion coefficient \xi) describes the relative motion 
between B (in state ^>i) and A. The intrinsic energy £j of B and the relative momen- 
tum Pi between B and A satisfy the energy conservation, Ei = h 2 Pf / (2/x) = E tot —ii, 
with the reduced mass fj, between B and A. 

Multiplying Eq. ([6]) by from the left, we can obtain a set of coupled differ- 
ential equations for \Xi), called CDCC equation, 

[Tr + tJi.i -E l ]\xi) = ~Yl Ui,i>\Xi>)- (8) 

The coupling potential E/iy is defined by 

IJu = {&i\U\£ v ). (9) 
The CDCC equation ([8]) is solved under the usual boundary condition for (R\xi) = 

2.2. Smoothing method 

In this subsection, we present a new way of smoothing discrete breakup S-matrix 
elements Si obtained by CDCC. The exact breakup T-matrix element to the exact 
continuum state ip(k) of B is given by 

T EX (k,P) = (i;(k),P\U\V) (10) 

with \ip(k), P) = \ip(k)) <g> \P), where k and P are momenta in the asymptotic region 
associated with the coordinates r and R, respectively; \ip(k)) is the exact two-body 
wave function of B with energy e satisfying 

[H B -sMk))=0, (11) 

and \P) is the plane wave function satisfying 

[T R - (E tot - e)]\P) = 0. (12) 

The exact three-body (A+b+c) wave function & can be replaced by the correspond- 
ing CDCC wave function i^ CDCC with good accuracy. Inserting the approximate 
complete set V of Eq. ([5]) between the bra-vector and the operator U of the right 
hand side of Eq. (|10p . we can obtain the approximate smooth T-matrix elements 
T(k,P): 

T(k,P) = ^mk)\£i)(£i,Pi\U\V c » cc ) 

i 

= J2^i(k)fi, (13) 
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where \<&i,Pi) = \<Pi) (g> \Pi), Fi(k) is the smoothing factor defined by 

Fi{k) = (^(fc)|*i>, (14) 
and Ti is the breakup T-matrix element of CDCC defined by 

t; = <^,^|t/|^ CDCC >. (15) 

Since the breakup T-matrix elements are proportional to the breakup S-matrix ele- 
ments, Eq. (fT3|) is reduced to 

S(k,P) = J2^i(k)Si. (16) 

i 

The smoothing factor J-i(k) is the overlap between the pseudostate <£j and the ex- 
act continuum state tp(k). We propose to use the complex-scaling method (CSM) 18 ) 
for the calculation of Ti{k). This new smoothing method is applicable to not only 
the two-body projectile but also the three-body one, as shown below. 

In CSM, the scaling transformation operator C(9) and its inverse operator are 
defined by 

<r|C(0)|/) = ^ 2e f(re ie ), (g\C~ l (9)\r) = {e^/ 2d /(re"")}* , (17) 

since 

{g\C-\e)C{6)\f) = J dr{g\C-\9)\r)(r\C{e)\f) = W) (18) 

for any normalizable states g and /, and hence C~ 1 (9)C(9) = 1. The transformed 
Scrhodinger equation is defined by 

(H^-s e )\^) = 0, $ e \(H$ i -e $ )=0 (19) 

with 

(r\H&\r') = {r\C(9)H B C-\9)\r') = 

where \x\, c is the reduced mass between b and c. The transformed Hamiltonian H° 
is not hermitian when 9 ^ and thus the ket vectors and the bra ones (ij) e \ are 
biorthogonal to each other. 22 ) Properties of and e e are summarized as follows. 
1. In the case that \ip e ) is a bound state \ifjj) with a negative energy e 6 - = 
—h 2 (kj) 2 /(2/j,b c ), (r\i/jj) has an asymptotic form exp(— k^re 10 ), so that it be- 
comes a normalizable L 2 -type function for < 9 < tt/2. Operating C(9) on 
Eq. (Illj) leads to (H^ — e)C(6)\ip) = 0. Comparing this equation with the first 
equation of Eq. (fT9j) . one can find that \ip®) = C(Q)\ij)j) and the eigenenergy 
is independent of 9. 



-2i6 



v r 

2fJ.bc 



+ V bc {re l 



8{r 



(20) 
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2. In general a resonance state \ip^) has an energy e e a = h 2 (k^e~ 1 ^) 2 /(2/j,f, c ) with 
positive 4> and kP a . The state has an asymptotic form exp(ik^re l( - 9 ~^), and 
hence it is a normalizable L 2 -type function for 8 > (p. Let us take the case that 
9\ > 6 > (j). Operating C{9\ — 9) on the first equation of Eq. (fT9j) and using 
C(0i) = C(9 1 -9)C{9), we obtain (H^ 1 -e e a )C{9 l -#)|<) = 0. Identifying this 
equation with (if* 1 - e e ^)\^) = leads to the fact that = C(9 1 - 6>)«) 
and e e a does not depend on 9 if > <fi. 

3. The scattering state \ijr{k)) should be normalizable as (tp 6 (k)\i{) e (fc')) = 5(k — 
k'). This means that the incident- wave part of the state is <po = (2ir)~ 3 ^ 2 exp(ik- 
r) and hence k has to be transformed by the scaling transformation as k — > 
ke" tS . Therefore, the scattering state has an energy e s (k) = (hke~ lS ) 2 /(2/i& c ), 
so that \^ e (k)) ^ C{9)\tp 0=o {k)). Note that (r\C {9)\ip 9=0 (k)) diverges at large 
r and thus not normalizable. 

The scattering state \ip{k)) at 9 = satisfies the integral equation 

\m) = l¥>o> + -^T^o) = |^o) + C- l (9)—^V b d c C(9)\ m ), (21) 
e — H-Q e - Hfc 

where use has been made of C~ 1 (9)C(9) = 1 in the second equality of Eq. (f2Tj) and 
\tpo) describes the plane wave satisfying 

[T r -e]\<po) = 0. (22) 

One can thus obtain the smoothing factor by 

= (## ) + (^IC-H^^t^^cC^I^o), (23) 

if the full propagator G® = (e — i^g)^ 1 is given for positive e. 

The eigenstates {tpj,ip^,,ip (k)} are normalizable and can form the complete 
set. 23 ) Thus, the spectral representation of G 9 is given by 

= I^I>7^<^I + ^ lO^^<^l+ / ^I^Cfc))^— (24) 

j j a a ^ ' 

In this representation, the resonance part (the second term of the right-hand side) is 
separated out from the scattering part (the third term) . This is useful in investigating 
a role of the resonance channel on the elastic scattering and the breakup reactions. 

The exact representation Eq. (|24p is still complicated and not so useful in partic- 
ular for systems more complicated than the two-body system. The full propagator 
(r\G e \r'}(= G e (r,r')) satisfies the equation 

G e (r,r') = Gl{ry) + f dr"G e (r,r")V(e id r")G e (r",r') 

= G e Q (r,r') + j dr"G e (r,r")V(e* e r")G e (r",r') + .-- , (25) 

where 

n 0, a / | 1 i a fi bc e 2ie exp[ik-\r-r'\e ie ] 

Go(r,r) = (r\ £ _ \r ) = j— . (26) 
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The free propagator Gg(r,r') is an exponentially damping function of \r — r'\ when 
< 9 < ir. Each term in the last form of Eq. (|25p has this property. Therefore, the 
full propagator G (r,r') keeps the damping property. This assures that G e can be 
expressed by a superposition of L 2 -type basis functions; namely, 

G e « V e G 9 V 9 , (27) 

where the operator V e represents the modelspace spanned by L 2 -type basis functions. 
The most convenient representation of V e is 

V° = (28) 

i 

where the <Pf are eigenstates of obtained by diagonalizing it with L 2 -type basis 
functions: 

(StlH^f,) = efSu, (29) 
One can find the approximate relation (^|G e |^,) ~ (e — ef)~ l 5n', because 

b lV = ffiffl) = {^\(e-H e B )G d \<Pl) « {^(e-H^G^) = (e-e? )(*i\G>\*fo. 

(30) 

Using this approximate relation, we can obtain 

g^ei^)(^i^x^i=Ei^)^(^i- ( 31 ) 

iV i 1 

This approximation has been applied to calculations of the continuum level den- 
sity, 24 ) the scattering amplitude, 25 ) and the strength function of the electromagnetic 
transition. 19 ^' 20 ) Validity of the approximate expression Eq. (|3ip in the CDCC frame- 
work will be justified, if the breakup S'-matrix elements converge as the modelspace 
V e is extended. Inserting Eq. (|3ip into Eq. (|23p . we obtain the equation for the 
smoothing factor: 

« (#^ o ) + ^^ T (^|C- 1 (0)|^)(^|yVf")(^l^)|^o). (32) 

i'i" £ £ i' 

When the Gaussian basis functions are taken, the matrix elements (<Pi\C~ l (9)\ <frj) 
and (d>f„\C (9)\ipo) are analytically obtained. The matrix elements {^,\V 9 \^ e i „) are 
also easily obtained by making a single integral. Thus, the final formula Eq. (|32p for 
the smoothing factor is quite useful. 

§3. Test calculations 

In this section, we test the validity of the smoothing formula, Eq. (|32p . i n two 
cases of the 58 Ni(c?,pra) reaction at 80 MeV and the 12 C( 6 He, 4 He 2 n) reaction at 229.8 
MeV. In the former, we have chosen d as a typical projectile with no resonance and 
treat it with the n+p model. In the latter, we have selected 6 He as a projectile 
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with resonance and for simplicity will treat it with the two-body 2 n+ 4 He model. In 
principle, 6 He should be treated as the three-body n+n-|- 4 He system. However, since 
our interest is the 2 + resonance in the present test, we use the dineutron ( 2 re+ 4 He) 
model that can describe the 2 + resonance reasonably well. 

The pseudostates <£j in the CDCC calculations are constructed by diagonalizing 
Hb with the complex-range Gaussian basis functions, which are found to effectively 
describe breakup processes. 10 ) The explicit forms of the basis functions are 

$(r) = /exp [-(r/a,) 2 } cos [vr/2(r/a,) 2 ] , (33) 

<fi e (r) = /exp [-(r/a,) 2 ] sin [7r/2(r/ aj ) 2 ] , (34) 
where j = 1—N. The range parameter aj is taken as a geometric progression, 

aj = a 1 (a N /a l )0'-i)/(^-i). (35) 

The number of the basis functions is 2N for each £. 

On the other hand, the complex scaled states 4>f are constructed by diagonalizing 
i?g with the conventional real-range Gaussian functions, 

</> jt (r) = r e exp[-(r/ aj ) 2 }, (j = 1-N), (36) 

which are valid for describing the approximate expression of G e as shown in Ref. 25). 

3.1. 58 Ni(d,pn) reaction at 80 MeV 

In the analysis of the 58 Ni(d,pn) reaction, we adopt the same model Hamiltonian 
as in the previous paper. 10 ) As for the breakup states of d we take i = and 2 
states and truncated them at k = 1.5 fin -1 . The range parameters of the complex- 
range Gaussian basis functions are (2iV = 40, a± = 1.0 fm, ajy = 30.0 fm). In the 
calculation of the complex scaled states, we take (N = 30, a,\ = 1.0 fm, = 30.0 fm) 
as the parameter set of the real-range Gaussian basis functions. 

Figure Q] shows the squared moduli of breakup S-matrix elements as a function 
of k at the grazing total angular momentum J = 17; panels (a), (b), (c) and (d) cor- 
respond to {I, L)=(0,17), (2,15), (2,17) and (2,19), respectively. Here L is the orbital 
angular momentum regarding R. The open circles are the results calculated with 
the exact smoothing factors that are obtained by solving the two-body Schrodinger 
equation (fTTj) numerically. The dotted, dashed, and solid lines correspond to re- 
sults of the new smoothing method with CSM at the scaling angle 6 = 5°, 10°, and 
15°, respectively. One can see that the results converge to the exact solution as the 
scaling angle 9 increases. 

3.2. ^CfHefHfPn) reaction at 229.8 MeV 

We take i = and 2 breakup states and truncate them at k =1.5 fm -1 . We use 
the following interaction V2 na between 2 n and 4 He: 

V 2na (r) = V exp[-(r/(3) 2 ], (37) 

where the depth Vq and the range are determined so as to reproduce the ground- 
state and resonance energies of 6 He. As for the interactions between 12 C and each 
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Fig. 1. The squared moduli of breakup S- matrix elements as a function of k at the grazing total 
angular momentum J = 17 for d+ 58 Ni scattering at 80 MeV. Panels (a), (b), (c), and (d) 
correspond to (£, L) = (0,17), (2,15), (2,17), and (2,19), respectively. 



constituent of 6 He, i.e. 6 He and 2 n, we take the optical potentials of d+ 12 C 26 ) 
and 4 He+ 12 C. 27 ) The parameters of the complex-range Gaussian basis functions are 
(2N = 40, a\ = 1.0 fm, oat = 30 fm), and those of the real-range Gaussian basis 
functions are {N = 30, a\ = 1.0 fm, = 30.0 fm). 

Figure [2] represents the squared moduli of the breakup S'-matrix elements at the 
grazing total angular momentum J = 25. Panels (a), (b), (c), and (d) correspond to 
(£,L) = (0,25), (2,23), (2,25), and (2,27), respectively. The peaks shown in panels 
(b), (c), and (d) represent the 2 + resonance of 6 He. The breakup 5-matrix elements 
calculated with the simple smoothing formula Eq. (|32p tend to the exact ones (the 
open circles) as 9 increases. One can conclude from Fig.[2]that the smoothing formula 
is valid also for the resonance part. 

§4. Summary 

In this Letter, we present a practical formula for smoothing discrete breakup 
5-matrix elements calculated by CDCC. This smoothing procedure based on CSM 
can easily be performed by diagonalizing H-q and the complex scaled one with 
the L 2 -type bases such as the Gaussian basis function and inserting the eigenenergies 
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Fig. 2. The same as in fig.[T]but for 6 He+ 12 C scattering at 229.8 MeV. The corresponding grazing 
angular momentum is 25. Panels (a), (b), (c), and (d) represent (£, L) = (0,25), (2,23), (2,25), 
and (2,27), respectively. 



and eigenstates into the formula. The validity of the formula is tested for two kinds 
of three-body breakup reactions, 58 Ni(g?, pri) at 80 MeV and 12 C( 6 He, 4 He 2 n) at 
229.8 MeV. For both the cases, the breakup S-matrix elements smoothed with the 
formula tend to the exact ones as the complex-scaling angle 9 increases. Thus, the 
formula is accurate and practical. In a forthcoming paper, we will investigate the 
practicability of this formula for four-body breakup processes of three-body projec- 
tiles such as 6 He and apply this method for analyzing the experimental data 28 ) -30 ) 
on the 6 He breakup reactions at lower and intermediate incident energies. 
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